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Spacecraft  Plasma  M&S 


*  Applications  from  Hall  thrusters  to  plumes  and  gas 
discharges 

*  Complex  physics:  excitation/ionization,  transport, 
radiative,  material,  etc. 


*  Multiple  spatial-temporal  and  density  scales. 


Current  focus: 

Develop  advanced  multiscale 
algorithms  for  plasma  M&S 
with  highly  non-equilibrium 
condition  and  collisional- 
radiative  kinetics 
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Collisional-Radiative  kinetics 


*  Non-equilibrium  modeiing  of  the  atomic  state  distribution 
function  (ASDF) 

—  Detailed  state-to-state  model  of  atomic  transition,  I.e., 
excitation  and  Ionization 

—  Rates  derived  based  on  ab  initio  cross  section. 

*  Exampies:  shocks 

*  Compiications: 

—  Numerical  stiffness  ,2 . . 


■  Atomic  CR 
Ar  shock  (31  levs) 
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Multiscale  problem 


Generic  Multiscale  problem:  ^  0 

Various  approaches  possible:  DNS,  multi-grid/AMR,  HMM, 
equation-free,  etc. 

—  Some  similarities  between  approaches 


Traditional  approach:  scale  separation  y  -  34  +  7/ 

—  when  £^0,  solve  ^  =  -fib's,  y/.x)  Flne-scale 

dX  £  4.  •  4. 

,  trajector- 

—  Use  "relaxed"  solution  for 

dx 

—  Works  slow  manifold  is  attractive  (neg.  eigenval 

Complications: 

—  Pos.  eigenvalues.  Ex:  instability,  inverse  casca 

—  Chaotic,  stochastic  fine  scales,  ... 


Coarse-scale 

trajectory 


—  Non-separation  of  scales 
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Multiscale  problem 


•  Generic  Multiscale  problem:  ■^  =  fo(y>x)+ e  ^  0 

•  Various  approaches  possible:  DNS,  multi-grid/AMR,  HMM, 
equation-free,  etc. 

—  Some  similarities  between  approaches 

•  Traditional  approach:  scale  separation  y  =  +  Xr 

—  when  £— >0,  solve  -fi(ys,yf,x) 

ax  E 

—  Use  "relaxed"  solution  for 

dx 

—  Works  slow  manifold  is  attractive  (neg.  eigenvalue) 


*  Complications: 


Atomic  CR: 

—  Pos.  eigenvalues.  Ex:  instability,  inverse  casc8^^  Hydrogen  (40  levs) 


—  Chaotic,  stochastic  fine  scales,  ... 

—  Non-separation  of  scales:  CR  kinetics 

15  orders  of 
magnitude 


10-^ 
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Kinetic  equation 


Boltzmann  equation: 

—  £  «1\  fluid  regime 

—  £~0(1):  kinetic  regime,  need  to  resolve  vfd 

dJ  +  v-V-J  =  -Q{f,f)+-. 

£ 


Q{f,f)=\\a{ 

d3 


V  -  V* 


>  (olf  ')  -  /(v)/(''*  y^codv' 


e(/,/)=o^/(v) 


—  For  collisional  plasma  (fully  and  partially  ionized):  need  to  include 
Coulomb  collisions  (FP),  excitation/ionization,  CX,  etc. 

Need  efficient  algorithms  for  collisional-radiative  processes  in  both 
regimes  and  transitional  regime. 
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Translational  Equilibrium 


*  Moment  method:  derived  by  taking  moment  of  the  kinetic  equations  with 
fiuid  ciosure 

—  5-moment  model  yields  Euler/NS  systems:  multi-species,  multi¬ 
temperature  CR  models. 

—  Magnetized  plasmas  are  often  modeled  with  MHD  with  a  hierarchy  of 
descriptions:  Ideal,  resistive.  Hall  MHD. 

—  Generalized  model:  multi-fluid  (Braginskii,  1965) 

*  CR  kinetics  =  rate  equations  for  each  excited  states 


—  Non-Boltzmann  population  of  ASDF 


—  Challenges:  many  states 


Atomic  hydrogen  (20  levs) 
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Translational  Non-equilibrium 


Not-too-far  from  equilibrium 

—  Discretized  EEDF  yields  rate  equations  for 
discrete  elements  (“bin”) 

—  DB  enforced  at  microscopic  level 

—  High-order,  implicit  and  energy  conserving 

—  More  efficient  compared  to  MCC. 

Far  from  equilibrium 

—  Particle  methods  with  DSMC  collisions 

—  Can  resolve  anisotropic  vdf 

—  Drawback:  slow  convergence,  reaction  branching, 
singular  rates,  particles  growth 


■excitation 


...  j  ... 

...  i  ...  bin  # 

Ui 

=  Ne  /  deUe) 

JEi 

dfii 

dt 

-  -Ni  fii  + 

j 

3 

dUj 

dt 

=  +iV,  5:  - 

i 

E%) 

i 
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Dynamical  Regime 


*  Examples:  thruster  plumes,  gas  discharges,  LPI,  etc. 

*  Require  hybridization:  fiuid/kinetic 

*  Key  chaiienges:  smooth,  efficient  and  consistent  transition 

*  What  we  need: 

—  Multiscale  statistics 

—  Coarse-graining  procedure  for  atomic  state 
—  Adaptive  fluid-kinetic  model 

*  Current  work: 

—  CR  level  grouping 
—  Particle  merging  schemes 
—  Multifluid  model 
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Level  grouping 


•  CR  modeling:  level-grouping  = 

•  Group  effective  rates  of  change  '  ' 

|j?i>n  ii^Fi  m  m<n  i^n  ”  j^m 

•  Internal  structure  of  group  is  assumed  Boltzmann  (7^) 

—  Piecewise  exponential 

•  This  does  NOT  mean  the  entire  ASDF  is  Boltzmann!! 

•  Group  temperature  must  be  determined  additional  conservation 
equation,  e.g.: 

•  Procedure?  * 


—  Solve: 


■a.jV.  [e  E 

LprL>n  i€n 

+  C.ipjT; 


jtm 


m<n  it™ 


with; 


^  =  W,  [AE) 


iterated 


computed 


C,{T„y.=  —{AE}^: 


dT^ 

¥ 


tabulated 


-However...  =  o(e) 

C,(T„)  ~  o(€) 


where  e  =  e 


6J:  = 
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Level  grouping 


CR  modeling:  level-grouping 

•  other  approaches? 

—  Sub-partitioning:  lowest  level  no  and  total  J%  (no  need  for  S^) 


J  V| 


=  io(c)/o(e)  ...fails 


—  Sub-partitioning:  lowest  level  no  and  upper  distribution 


rp*2 

~ - -a — 

"  -  Zi,{T*){AE)ATi 


■9^ 


=  o(t) /£?(£}  ...fails 


—  Approximate  Z^,  by  expanding  around  mean  energy:  ~ke„  = 

^ ^  |l  +  . . .j  where  6i  =  AE, --KE„ 

•  With  no,  partitioning:  1/ln(1+s)  ...fails 

•  With  no,  partitioning:  1/ln(s)  ...succeeds! 

•  Improve  with  successive  iterations...  .  _  -p  k-t  'i 
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Level  grouping:  OD  test 


*  Isothermal  heat  bath 


IQ-U  . I 

0  2  4  6  8  10  12  14 

£  (eV) 


Works  very  well  (also  tested  in  cooling  regime). 
Much  better  than  uniform  (standard)  grouping. 
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Level  grouping:  Energy  conservation 


*  CR  modeling:  level-grouping 

•  Look  at  uniform  bins  -  2  formulations: 


jera  "  Energy  conserving! 

ji=m  ieri 


Right-form  is  product  of  Energy  x  dJ7Vn/dt... 

—  Conservation  follows  from  definition  £„  = 

•  Same  principle  can  be  applied  to  Boltzmann  groups: 
—  Start  with: 


=  {A£) 


dt 


itn 


—  Express  in  terms  of  conserved  variables: 


{fTn 

IT 


dt 


—  ulnr 


dt 


with 


= 


a.nrX? 


J  _  g,  A'nr 

and 


no 
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Level  grouping:  Energy  conservation 


*  CR  modeling:  level-grouping 

•  Finally... Procedure  shown  to  be  equivalent  to  replacing 
energies  by  “effective”  (condition-dependent)  values  («  EOS) 

Excitation: 


|nf )  ~  ~  j  hO 


with 


Enf  =  £’^-1-  +in> 

10-11 


NOW,  energy  is  conserved 

(down  to  round-off) 


10 


-II] 


10 


-la 


l-l 


10 


-M 


10 


-IJS 


H]r^> 


icr®  icr’' 

time  {see) 
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Particle  Merging  Schemes 


The  schemes  consists  of  merging  3 
or  more  particies  to  2.  Mass, 
momentum  and  kinetic  energy  are 
exactiy  conserved;  Electrostatic 
energy  also  conserved  in  physicai 
space. 

Split  analogousiy  defined  by  merging 
oniy  fractions  of  original  particles. 

To  inhibit  thermaiization,  an  octree  in 
velocity  space  is  used  so  that  oniy 
near  neighbor  particies  are  merged. 

Higher-moment  conserving  schemes 
obtained  with  increased  number  of 
merge  resuit  particies  generated. 
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Examples 


•  from  OD  ... 
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Timej  (ns] 


Examples 


•  from  OD  to  1D3V... 


—  Gas  discharge  with  ionization. 


Density  (m‘^) 


l,«+0®  1^+10  l,e+ia  l.e+14 


Computational 
Particles  per  Cell 


1.  l.e-l-Ol  1.e+02  l.e+D^  l,c-l-04 
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Examples 


•  from  OD  to  1D3V  and  3D3V  ... 

—  Annular  test  case 
—  More  test  cases  underway 
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Normalized  velocitj^ 


Multi-fluid  equations 


*  Multi-fluid  model: 

—  Develop  series  of  models  for  variable  conditions: 

Electrostatic  +  Magnetostatic  ^^Electro-magnetic  3-fluid 


—  Use  implicit  models  to  eliminate  constraint  of  sequence 

of  fast  time  scales:  c,  v^,  ^ 

—  Price  to  pay:  lack  of  resolution  misses  physics... 


Theory  -  =  22.7  rad /nsec 

Elxplicit  — I — 

ImpUcit  — • — 


0  0.5  1  1.5  2  2.5  3 


Transit 
Ionization 


Ion  Cyclotron 


10" 


Ion  Plasma 
Electron  Total 


10® 


Ion  Neutrsil 


Electron  Plasma 


10'"  10"  10® 
I 

Lower  Hybrid 


-4-  I  '> 


Excitation 

Electron  Ion 


10' Hz 
I 

Electron 

Cyclotron 


time  (nsec) 
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Multi-fluid  equations 


*  Multi-fluid  model: 

•  5-moment: 


dtPs  -\-  V  •  (psUs)  — 


dtiPs^s)  +  V  ■  (fsUsUs  -I-  Ps)  =  ZsCRs  (E  -I-  Us  X  B)  -|-  f " 


9t£s  "H  V  ■  {Us^s  -\-  Ps  ■  Us)  +  V  ■  C|s  — js  ■  E  +  Wj 


•  Add  Maxwell’s  equations:  „  ^  ac  e 

V  X  E  =  -^tB  V  -  E  =  — (ZiHi  -  n.) 

V  xB  =  ^Qj-  -^9t£  V  .  B  ^  0 

•  Add  collisions: 

—  Elastic  -  Bragiinski  terms 

—  Inelastic  -  warning!  Rates  depend  on  both  T  and  relative  velocity 


fnfegC^US'^  ^i,^2)d^ldn2d^Vnd^V, 


—  Multi-fluid  CR  model  from  fundamental  principles  being  developed 
(incl.  detailed  balance) 
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Multifluid  equations 


*  Electromagnetic  shock:  generalized  Brio-Wu^ 

—  FV  with  WENO  reconstruction  and  RK3 


^Shumlak  &  Loverich,  JCP  2003 
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Laser-Plasma  Interaction 


*  Ion  acceleration  due  to  ponderomotive  force 


X  |/im| 
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Conclusion 


*  Multiscale  algorithms  for  nonequilibrium  flows  with  CR 
kinetics 

—  Particle  merge/split  for  particle  management,  efficient 
sampling,  inelastic  collisions  ... 

—  Level  grouping  schemes  of  electronic  states,  for  dynamical 
coarse-graining  of  ASDF. 

—  Multifluid  equations  to  efficiently  capture  electron 
“hydrodynamics” 

*  Ongoing  works: 

—  High-order  particle  merging  schemes 

—  Multi-D  simulation  with  level  grouping 

—  Modeling  of  inelastic  collisions  in  multifluid 
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